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^ , Due to the strongly nonlocal nature of fxc{r,r',uj) the scalar exchange and correlation 

QP ' (xc) kernel of the time-dependent density-functional theory (TDDFT), the formula for 

' Q the friction coefficient of an interacting electron gas (EG) for ions tends to give a too 

large value of Q for heavy ions in the medium- and low-density EG, if we adopt the local- 
density approximation (LDA) to /a;c(r, r', oj), even though the formula itself is formally 
■ exact. We have rectified this unfavorable feature by reformulating the formula for Q in 

' terms of the tensorial xc kernel of the time dependent current-density functional theory, 

' to which the LDA can be applied without intrinsic difficulty. Our numerical results find 

themselves in a considerably better agreement with the experimental stopping power of 
Al and Au for slow ions than those previously obtained within the LDA to the TDDFT. 



1. Introduction 

The problem of the stopping power (SP) of solids for moving ions (the energy-loss 
by ion per its unit path-length) has been attracting much attention for decades (see, 
e.g., Ref . for recent reviews). Theoretically, this problem has been investigated 
in two distinct ways: One is a perturbation approach in which the projectile-target 
interaction is taken as an expansion parameter^. The other is based on the density- 
functional theory (DFT), treating the problem in a fully nonlinear manner from the 
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outset sua. 

For slow projectiles, the perturbation expansion generally fails making the 
nonlinear treatment indispensable. In the limit of zero velocity, the potential scat- 
tering (binary-collision) theory is known to give an expression for the friction coef- 
ficient (FC) (SP divided by the projectile velocity at its zero value) of electron gas 
(EG) as El IS HI 

Q ^ nokFatrikp), (1) 

where kp is the Fermi momentum, atrikp) is the transport cross-section of the 
elastic scattering of an electron at the Fermi level in the self-consistent Kohn-Sham 
(KS) ^ potential VKs{f) of the ion statically immersed in the EG, and fiQ is the 
density of the homogeneous EG in the absence of the ion. Equation ([1]) is, however, 
incomplete El El in the sense that it ignores the role of the dynamical exchange and 
correlations (xc) effects in the problem of the ion slowing. These effects cannot be 
neglected even at vanishingly small velocities. The complete formal solution to the 
problem of the friction coefficient of the interacting EG including the dynamical xc 
effects has been recently found as 

Q = Qi + Q2, (2) 



where 



V 

dlmxKs{r,r',uj) 



=0 



dbJ 

[Vrno(r) ■ v] [Vr'no(r') • v] 



drdr', (3) 



9Im/j;c(r,r', w) 



. dvdv'. (4) 

In Eqs. ([2])- 11]), Qi and Q2 denote, respectively, the independent-electrons and the 
dynamical xc contributions to the friction coefficient; x/f5(r, r', w) and /a,c(r,r',w) 
are the KS density-response function and the dynamical xc kernel ^1^, respectively, 
of the inhomogeneous many-body system of an ion at rest in the EG, no(r) is the 
ground-state density of this system, and v is a small velocity of the projectile, on 
the value of which the final results evidently do not depend. 

The dynamical xc kernel is defined as the Fourier transform with respect to time 
of the functional derivative of the time-dependent xc potential with respect to the 
time-dependent electron density 

^T4cW(r, t) 



/,eK(r)](r,i;r',i') 



(5) 

n=no(r) 



5n{v',t') 

A formal proof of Eqs. ©-(HI) is provided in Appendix [Bl 

For the independent-particle contribution, Eq. ^ can be shown to coincide 
exactly with Eq. ^ of the binary-collision theory (see Appendix |D] of the present 
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paper). Actual evaluation of Eq. ([T]) can be done through a well known procedure 
of the self-consistent solution of the KS equations 

[-(1/2) A + VKsir)]iJ^{r) ^ e,ij,{r), 

where ipiir) and are the single-particle wave-functions and the eigenenergies, 
respectively, in the potential 

Vxsir) = -Zi/r + J [no(r') - no]/\r ~ r'| dr' (6) 
+ yxc{nQ(T)\ - 14c (no), 

where Z\ is the atomic number of the ion, and V^c[no(?')] is the xc potential of the 
static DFT. The ground-state density is determined as 

no(r) = I^^WI'' (7) 

where e_p is the Fermi energy. Upon achievemen t of self-consistency, the transport 
cross-section can be conveniently found through 

AfT 

;=o 

where 5i{k) is the phase-shift of the scattering state with momentum k and angular 
momentum /. 

Determination of both the noninteracting-electron (Qi) and the dynamical xc 
{Q2) parts of the friction coefficient require, respectively, pertinent approximations 
to the static xc potential Vxcir) and the dynamical xc kernel /xc(r, r', w). The 
local-d ensity approximation (LDA) provides simple and rather accurate schemes 
[131 ttH for the quantities derived from the static xc potential. The situation with 
the dynamical xc kernel is, however, much more complicated: The scalar xc kernel 
/a;c(r, r', a;) of the time-dependent density functional theory (TDDFT) is known 
to be a strongly nonlocal function of space coordinates 1^^, implying a difficulty 
in directly constructing an accurate approximation to calculate Q2 through Eq. 
(j4]). On the other hand, it is known that if the current-density rather than the 
particle-density is chosen as the basic variable of the theory, then the corresponding 
tensorial xc kernel is a function of short range. This fact made it possible to construct 
a consistent LDA sche me with in the framework of the time-dependent current- 
density functional theory tEl Ull (TDCDFT). More recently, the TDCDFT has been 
successfully applied in studying the damping of excitations in disordered systems 
and dynamical corrections to the resistance of point contacts 
Therefore the main purpose of this paper is to recast Eq. ^ into a form ex- 
pressed through the tensorial xc kernel of the TDCDFT which allows the LDA 
treatment. With using this form, we have made comparison between TDDFT and 
TDCDFT through numerical investigation to find that Q2 is overestimated in the 
LDA to the TDDFT for heavy ions immersed in the medium- and low-density EG, 
while the LDA to the TDCDFT rectifies this unfavorable feature. 
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It must be noted, that the exact TDDFT and the exact TDCDFT, were they 
possible to develop, would yield, of course, exactly the same friction coefficient Q. 
Both theories would then involve nonlocal scalar and tensorial xc kernels, respec- 
tively. The basic idea of our method is to express quite generally the scalar xc kernel 
of TDDFT through the tensorial one of TDCDFT [Eq. ^]. Then we show that 
if the local tensorial xc kernel of TDCDFT is used on the right-hand side of this 
expression, the resulting nonlocal scalar xc kernel of TDDFT is free of the inconsis- 
tencies of the LDA to TDDFT, and when used in Eq (jl]), it yields a good description 
of the many-body dynamical effects in the friction coefficient of the interacting EG 
for ions. 

In Refs. Uni EOl j^g^^g summarized main resuUs of the TDDFT and TDCDFT 
approaches, respectively, to the problem of the stopping power of metals f or s l ow 
ions. The present work gives rigorous proofs of all the assertions made in Ref. 
as well as specifies details of both the mathematical formalisms and the computa- 
tional procedures. 

The organization of this paper is as follows: In Sec. [2] and Appendices A - ID] 
we present a formal TDDFT of the friction of ions in EG and in Subsection 12. II we 
discuss a fundamental difficulty the LDA to the dynamical xc kernel encounters, 
i.e., the violation within LDA of the requirement that the friction coefficient of free 
space for isolated atoms and ions must be zero. In Sec. [3] we develop the TDCDFT 
of the friction of ions in EG and demonstrate the resolution within the LDA to the 
TDCDFT of the above contradiction (Subsection 13 . 21 and Appendix IE)) . In Sec.[4]we 
outline the calculational procedures, present the results and their discussion. Our 
conclusions are collected in Sec. [5] Appendix [F] provides explicit expressions for the 
response functions and xc kernel in the spherical coordinate system that we use in 
our calculations. 



2. Time-dependent density-functional theory of the friction 
coefficient 

We consider a recoiless probe particle of the charge Zie, where e is the absolute value 
of the charge of electron, moving with velocity v in the otherwise homogeneous gas 
of interacting electrons at zero temperature. The stopping power is the retarding 
force that the polarization charge distribution in the vicinity of the projectile exerts 
on the projectile itself. Accordingly, one can write^^ (we use atomic units throughout 
the paper) 



dE Z f 

--T = / drdr'(5(r-vt)v-Vrn„d(r',i)/|r-r'|, 

ax V J 



(9) 
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Uind being the electron density induced by the projectile. Our first step is to obtain 
the following expression for the friction coefficient of the interacting EG: 

g^lini-i^^ (10) 
v^o V ax 



[VrK)(r) ■ v] [Vr'Voir') ■ v] 5Inix(r, r', c^) 



dr dv' , 

j=0 



V V duj 

where Va{r) = —Zi/r is the bare Coulomb potential and xi^j^'i^) is the linear 
density-response function of the system of interacting electrons in the static field 
of an impurity of charge Zi at r = 0. A proof of Eq. PH)) is provided in Appendix 
A. We emphasize that Eq. (|10p is fully nonlinear with respect to the interaction of 
the bare charge Zi with the EG. The next step, which is performed in Appendix 
m is to show that Eq. ^ leads to Eqs. © - ©• The equivalence between Eq. (O 
and the friction coefficient of Eq. ([T|) from the binary-collision theory is proven in 
Appendix |Dj 



2.1. LDA to the TDDFT and the violation of the requirement for 
the friction coefficient of free space to be zero 

Within the LDA, the scalar dynamical xc kernel is given by El 

/,e(r, r', oj) = /,^, Jno(r), c.] 5{v - r'), (11) 

where /^^ jin^iS) denotes the q limit of the frequency-dependent longitudinal 
xc kernel of a homogeneous EG of the density n. Due to the spherical symmetry of 
no(r), substitution of Eq. pT|) into Eq. ([4]) yields 

n ^vr 7 ,i2 ^Im/,\^K(r),o;] 

(32 = — ^ / c^'^r^W 75 • 12 

6 J auj w=o 



An important test of the theory is the limit of zero density of the EG, i.e., 
no — *■ 0. It is evident that the friction coefficient for any atom (ion) should be zero 
in this case. The noninteracting-electrons part Qi of Eq. ^ obviously satisfies this 
requirement. However, it is easy to see from Eq. (|12p that within the LDA, Q2 
remains finite since the ground-state density no(r) of an isolated ion (atom) has a 
nonzero gradient, and the frequency derivative of Im/^^^(n,a;) is always nonzero 
and negative Therefore, the LDA to the scalar xc kernel does not satisfy the 
requirement of the friction coefficient of vacuum for an atom to be zero. 

Let us now demonstrate that the rigorous formula (|4]) does pass the above test. 
For a n ar bitrary bounded (finite) system, which is an atom in vacuum, a sum rule 
holds HSl 

j f,c{r, r', Lu) Vr'uoir') dr' = V,Kc(r). (13) 

Substituting Eq. (fT^ into Eq. ^ and noting that the right-hand side of the former 
does not depend on frequency, we obtain zero identically, which proves our assertion. 



20, 2008 2:48 WSPC/INSTRUCTION FILE arXiv 



6 

We note that the sum rule (|13p does not hold for extended systems and, therefore, 
Q2 is finite for atoms (ions) in EG of nonzero density. 



Table 1. Inaccuracy of LDA to TDDFT: Friction coefficient of 
free space (r^ = 00) and that of an electron liquid of = 2.2 for 
several atoms. Line 3 is the ratio of lines 1 and 2 (%). 



Atom 


He 


Be 


C 





Ne 


Mg 


Si 


Q{rs = 00) 


0.04 


0.11 


0.17 


0.24 


0.30 


0.36 


0.43 


Q(rs=2.2) 


0.34 


0.43 


0.70 


0.46 


0.16 


0.15 


0.54 


% 


12 


25 


24 


52 


188 


240 


80 



To give an idea how large can be the error introduced by LDA to TDDFT, in 
Table 1 we list the values of EC of free space (r^ — oo) for a number of atoms. For 
comparison, in column 3 of this table we list the EC of EG of rg —2.2 for the same 
atoms obtained with use of Eq. ^ (neglect of dynamical xc). 

The shortcoming of the LDA to the TDDFT manifested in the the violation of 
the requirement of zero friction coefficient of vacuum will be overcome in the next 
section within the LDA to the TDCDFT. 

3. Time dependent current-density functional theory of the 
friction coefficient 

The purpose of this section is to express Eq. ([4]) for the dynamical xc contribution to 
the friction coefficient in terms of the xc tensor of the time-dependent TDCDFT in 
order to overcome the difhculties the LDA encounters when app lied to the ordinary 
TDDFT as discussed in the previous section. We know that LtSi 

/,,,,,(r,r',^) = Xifs,.j(r,r',^) - X^\r,r',u;) + -^V,^^^ V^, (14) 

where fxdj is the tensorial xc kernel of the TDCDFT, Xij is the many-body cur- 
rent density-vector potential response function, and XKS,ij is the single-particle KS 
counterpart of the latter. On the other hand, the scalar xc kernel of the ordinary 
TDDFT can be written as 

fxc{r,r',uj) = Xxs('"''"''^)-X"\i",i"',tj)- 1^^^-^, (15) 

where x is the longitudinal density-scalar potential response function and xks is 
its Kohn-Sham counterpart. We can write in operator notations 

X = -4v-x-V, (16) 

and consequently 

X-^ = -^V-2V.(ix^)"'-VV-2, (17) 
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where L is the longitudinal projector operator 



Using a simple operator identity 



(18) 



where T = 1 — i is the transverse projector, we can write for the inverse scalar 
response function 

, ,2 



X 



X. 



-'^x^fr'Tfr' 



•VV" 



(19) 



Then by Eq. we have 



c 



{/.c+( 

(xK\-hc)f\ '\xK's~hc) 'Xk's {txksT) "'xKsj • VV 



(20) 



Equation ((20|) constitutes a formal expression for the scalar xc kernel fxc of 
the ordinary TDDFT in terms of the tensorial xc kernel fxc of the TDCDFT and 
the independent-particle (KS) tensorial response function xks- Then by virtue of 
Eqs. (HOI), dH), (Ull), and ^ one can write 



Q2 



^ cl"?o ^ - nQ)v-Llfxc+ [xks - /^c j r {xks ^ f^c) ^ 



{xks - f^<) - Xks (tXksT) Xks^ ^("o - "o)v. 



In the above equation, the operator in the braces is longitudinal from both sides, the 
fact of which can be verified by applying the T operator from cither side yielding 
zero identically. We therefore can drop L operators to the result 



Q2 



- lim uj Im (no - no)v • I f^c + [xks - fxc 



T 



(xks f^ 



T 



(21) 



Equation (|2ip is a formal expression for the dynamical part of the friction co- 
efficient in terms of the TDCDFT quantities. It, however, is not easy to use in 
calculations. We, therefore, proceed to derive an equation which will be shown in 
Sec. m convenient to implement in practical calculations. First, we introduce the 
notations 

X - (txk'sT) Xk's ■ i^o - no)v, (22) 

(23) 



y={Tx-'T) r •K-no)v. 
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Then this is straightforward to show that y satisfies the equation 



y = ^+(fxK'sT) /.c • [y - K - "o)v] 



Finally, with use of Eq. ([24| , Eq. (j2ip can be rewritten as 

Q2 = ^ lim uilm [x - {uq - no)v] ■ f.^c ■ [y - {no - no)v]. 



(24) 



(25) 



3.1. The Local Density Approximation for the tensorial xc kernel 

Within the LDA to the TDCDFT one can write 

fxc.iji'r^v' ,uj) Sj{v')dr' = — x 

rALDA 



V,y,^,^^^(r, iv) + — -V, a,,,.,(r, u:) 
no(rj 



where 



1 



y/,^^-^(r,c.) = -6;'Jno(r)]V,s,(r), 



lUJ 

^xc{n) is the xc energy density of the homogeneous EG of density n, 
o-a;c,ji(r,w) = ryj;c[no(r), [VjUi(r) + ViUj(r) 
2 , 



■ VfcWfc(r) 4 



+ Cxc[no{v),uj\ VkUkir) 5i 



where 

u(r) = s(r)/no(r), 
the viscosity coefficients are given by 



(26) 



(27) 



(28) 



(29) 



and /^^ ^(rt, ui) and /^^ j.(rt, w) are, respectively, the longitudinal and transverse xc 
kernels of the homogeneous EG with the density n. 



3.2. Resolution within the TDCDFT of the problem of the finite 
friction coefficient of vacuum 

To show that within the LDA to the TDCDFT the friction coefficient of vacuum 
is zero, it is sufficient to prove that with the tensorial xc kernel of the previous 
subsection and the scalar xc kernel obtained from it by Eq. ([^0]) . the sum rule of 
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Eg . p3|) holds. We proceed by recalling exact sum rules for the tensorial quantities 



fxc.ij{r,r' ,uj)no{r')dr' = 7 ViVjV^c(r), 



(30) 



J x^kir,r',uJ)V',V^Voir')dr' J (r, r', c^) dr' - no(r) <5y , (32) 

where Vxc{j^) is the static xc potential. With the LDA to fae of the previous section, 
the sum rules ([5D|) - (15^ are satisfied by construction In Appendix |E] we prove 
that for a finite system the tensorial sum rules ([30|) - ([32|) entail the scalar sum rule 
(fT3|) . As soon as the latter is satisfied, results of Sec. 12.11 lead to the zero friction 
coefficient of vacuum within the LDA to the TDCDFT. 



4. Calculational procedures, results, and discussion 
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03 
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■q 
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o 
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2 4 6 8 10 12 14 16 18 20 



1 

Fig. 1. Friction coefficient of a homogeneous EG of r'a = 2.2 versus the projectile charge Zi. 
The open squares are the results of the calculation with neglect of the dynamical xc as obtained 
from Eq. JTJ. Triangles are the results of the calculation with the dynamical xc included within 
the LDA to the conventional TDDFT as obtained from Eq. 11121 . Open circles are the results of 
the calculation with the dynamical xc included within the LDA to the TDCDFT as described in 
Sec. |4] Solid squares are the measured SP of Al of Ref. ESlfor slow ions (v = 0.5 a.u.) moving at 
a distance of 1.2 a.u. from the last atomic plane of the Al (111) surface. 
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Fig. 2. Friction coefficient of a homogeneous EG of = 2 versus the projectile charge Zi. Open 
squares are the results of the calculation with neglect of the dynamical xc as obtained from Eq. |[T||. 
Triangles are the results of the calculation with the dynamical xc included within the LDA to the 
conventional TDDFT as obtained from Eq. I|12| l. Open circles are the results of the calculation 
with the dynamical xc included within the LDA to the TDCDFT as described in Sec. |4] Solid 
squares are the measurements of Rof. of the SP of Au for slow ions (v = 0.68 a.u.) channeled 
along the (110) direction. The dotted line is the calculation of Ref. El with the dynamical xc effects 
included in the framework of the linear-response theory. 



We solve Eq. p4)l for the y vector at finite frequencies, then substitute the 
resuhs into Eq. ((25)l . The friction coefficient is found by the extrapolation from 

finite to zero frequency. To conveniently treat ^TxJ^^gT^ in the right-hand side 
of Eq. (124]) , we use the identity 

(^X/fS^) = XA'S + ^XKS ■ Vxk5^ ■ XKS, (33) 

and equations (|F.2P - (jF.5P and (jF.7P - (jF.8P are utilized for the explicit evaluation 
of the xks and f^c operators, respectively. With these provisions, the numerical 
procedure is to approximate Eq. (|24p by using a complete ortho-normal set of basis 
functions thus reducing the problem to a system of linear equations. The natural 
choice for the radial basis functions is 

Mr) = [a^{n + l)in + 2)]~'^^ I L^^\r j a) , (34) 
where L^n^ {x) are the associated Laguerre polynomials 0, and a is a scaling param- 



'^k = 2 ensures the orthonormality of the radial basis functions with the weight of . 
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Fig. 3. Friction coefficient of a homogeneous EG of carbon density (rs = 1.59) versus the pro- 
jectile charge Z\ . Open squares are the results of the calculation with neglect of the dynamical xc 
as obtained from Eq. Q. Open triangles are the results of the calculation with the dynamical xc 
included within the LDA to the conventional TDDFT as obtained from Eq. II12I I. Open circles are 
the results of the calculation with the dynamical xc included within the LDA to the TDCDFT 
as de scribed in Se c. Bl Solid circles, triangles, and diamonds are transmission measurements from 



Refs. IIH and IMof the random SP of C for ions moving with velocities v = 0.41, 0.83, and 
0.25 a.u., respectively. 



eter. Since the exact results should not depend on the choice of the basis functions, 
we have been checking the convergence of the numerical procedure by varying a 
and making sure that the respective results differ insignificantly. 

In Fig. [Tl we plot the friction coefficient of EG of = 2.2 versus the atomic 
number of a moving ion. Results of the calculations with neglect of the dynamical 
xc [the binary-collisions approximation, Eq. the LDA to TDDFT [the sum of 
Eq. HI) and Eq. and the LDA to TDCDFT [the sum of Eq. ^ and Eq. (|25| ] 

are shown together with the experimental data of Ref. for ions moving with 
the velocity of 0.5 a.u. at the distance of 1.2 a.u. from the last atomic plane of the 
(111) surface of aluminum. The inhomogeneity of the electron density the ions travel 
through is not strong under these conditions, and in the calculations we have used 
Vs estimated experimentally Moreover, the experimental SP is predominantly 
electronic since the trajectory of the projectile ion remains well separated from the 
lattice atoms. Together, these two conditions justify the comparison with the theory 
within the homogeneous EG model. The non-monotonic dependence of the friction 
coefficient on the atomic number of the projectile (the so-called Zi-oscillations) is 
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known to result within the single-particle theory from the competition between the 
increase in the EG-ion interaction with the growing charge of the bare nucleus of the 
ion and the decrease of the same interaction due to the screening by the formation 
of shells of bounded electrons of the pseudo-atom as well as its resonant states, as 
discussed in detail in Ref. SI. 

While the LDA to the conventional TDDFT (triangles in Fig. [T]) largely over- 
estimates the friction coefficient at Z\ > 5, the LDA to the TDCDFT (open cir- 
cles in Fig. [1]) is in a good agreement with the experiment in a wide range of 
3 < Zi < 18. The significant deviation of our results from the experiment occurs 
at Zi = 1,2,19,20, where the experimental friction coefficient is lower than the 
binary-collisions calculations (open circles in Fig.[T]). This feature has been recently 
shown to be a single-particle effect due to the finite velocity of the projectiles 
hence it is an effect of the deviation from the linear dependence of the SP on the 
velocity. The same effect gives positive contribution at 8 < Zi < 12 suggesting that 
combined with the many-body effects of the present theory the agreement with the 
experiment can be further improved. El In the range 13 < Zi < 17, the dynamical 
many-body effects seem to be solely responsible for the enhancement of the friction 
coefficient compared with the binary-collisions calculation results. 

In Fig. [21 we plot the friction coefficient of EG of Ts = 2 versus the atomic number 
of ions in the range 5 < Zi < 39. This is compared with the available measured SP 
for ions with the velocity of 0.68 a.u. channeled along the (110) direction in gold. 
Now because of the channeling, the collisions with the lattice atoms again do not 
give significant contribution to the SP. We can see that the general trend is that the 
LDA to TDCDFT improves the agreement between the theory and the experiment. 
However, within the range 16 < Zi < 19 the dynamical xc contribution is too 
small to account for the onset at the experimental data, nor can the persistent 
enhancement of the friction coefficient in this r ang e be attributed to the finite 

indi 

velocity effects within the single-particles theory Further studies are required 
to elucidate the nature of this onset, one of the possible causes being evidently the 
band-structure effects. 

In Fig. [21 we plot the results for EG of Ts = 1.59 corresponding to the valence 
electron-density of carbon. It is instructive to note that within the range 1 < Zi < 14 
the resuhs within the LDA to the TDCDFT and the LDA to the TDDFT do not 
differ significantly, which can be believed to be true in the general case for light 
atoms in the high-density EG. Then at higher Zi rather abruptly the dynamical 



''There is no reason, of course, for these two effe cts t o be additive, ft must be noted, that Ref. 
incorrectly attributes the overestimation by Ref. ESI of the contribution of the many-body effects 
to the use of the total ground-state density rather than the density of the delocalized states only. 
The total particle-density is, however, the basic variable of the TDDFT and without any artificial 
assumptions it enters the ri goro us result of Eg. l|4]l. The real source of the overestimation of the 
dynamical xc effects in Ref. was, as Ref. had anticipated and the present work shows, use 
of the LDA within the conventional TDDFT. The present work is overcoming this shortcoming 
within the LDA to the TDCDFT. 
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xc contribution almost vanishes, which effect also can be understood qualitatively 
recalling that for heavy atoms immersed in EG the charge-density distribution (the 
main quantity entering the theory within the LDA to the TDCDFT) is close to that 
of isolated atoms, and hence according to the results of Sec. l3.2l it should disappear. 
The experimental data in Fig. [3] correspond to the SP for ions moving along random 
trajectories inside bulk carbon and, therefore, they are strongly influenced by the 
collisions with lattice atoms, which fact hinders the quantitative comparison with 
the calculations within the EG model. 

It must be noted that the LDA to TDCDFT has no strict iustification in the 
Lo ^ limit.^^ Nevertheless, the calculations of the static polarizabilities of con- 
jugated polymers using the Vignale-Kohn functional in the form of Eqs. ([27| - (|29p 
have proved to be very successful .1221 EHI this respect, the situation seems to be 
analogous to that with the LDA to static DFT, which had been justified for weakly 
varying electron densities only, while the range of its successful applications has 

on on 

proved to be incomparably broader. Similar to the method of Refs. ^ we have 
been performing our calculations at finite to then extrapolating to cj = 0. 

5. Conclusions 

Within the current-density functional theory, we have derived equations for the 
dynamical contribution to the friction coefficient for ions moving in electron gas. The 
resulting theory demonstrates a crucial advantage over the ordinary time-dependent 
density-functional approach since within the local-density approximation the former 
ensures the correct zero friction coefficient in the limit of zero density of the electron 
gas, while the latter does not. 

We have traced the breakdown of the LDA to the ordinary TDDFT to the 
violation of a previously known sum rule for the scalar exchange and correlation 
kernel. 

Finally, we have performed calculations within the framework of the new theory 
with the results finding themselves in a good overall agreement with the exper- 
imental electronic stopping power of (i) the near-surface region of aluminum for 
slow ions moving in the grazing geometry and (ii) gold for slow ions moving under 
a channeling condition. 

Appendix A. Proof of Eq. (fTOj) 

For a projectile moving in a homogeneous medium, we can rewrite Eq. ^ in the 
reciprocal space as 




where q — (q, u) and 



27r2q2 ■ 
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An external potential 4>ext{q) induces in the homogeneous EG the particle density 

oo „ 

riindiq) = X] / d-qi-dqn-iXn{q,qi, ■■■,qn-i)<l>ext{q - qi) 



n=l ■ 



X<t>ext{qi - q2)--4>ext{qn-2 - qn-l)<t>exi{qn-l), 

where Xn{q,qi, ■■■qn-i) is the n-fold many-body density-response function of the 
homogeneous EG which satisfies n — 1 symmetry relations 

Xn{qo, ■■■qj,qj+i, ■■■,qn-i) = 

Xn{qo,---qj,qj-qn-i:---,qj-q:f^i), < j < n - 1. (A.l) 
Since the bare potential by the projectile is 

<Pext{q) =Vgd{uj - V -q), 

we can write 

oo „ 

Tiindiq) = S{ui -V ■ q)'^ / Xn{q,qi,-,qn-i) 

n=l •' 

where now q = (q, v • q). Then we can write for the stopping power 

oo 

V • qo 



xfn{q,qi, ...,qn-i) dq_o---dq_n-i, (A.2) 
where we have introduced the notation 

fn{qo,qi,-,qn-i) = Xn{qo,qi,-,qn-i) 

For the friction coefficient, which is the stopping power divided by the projectile 
velocity at its zero value, we can write 

oo n— 1 



n=l 7=0 



X /nj(qo,-,qn-i)rfqo-"rfqn-i, (A.4) 

where we have introduced further notations 

dfn{qo,-,qn-i) 



Our purpose is to prove that 



a;o=0,.-MWn-i— 



Q = P, (A.5) 
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where 



(27r)3lm / dqdk^^^^^^v.Vk 

J V V ^ 

ax(q,k,w) 



(A.6) 



du) 

and x(q, k, uj) is the Unear density-response function of the system of the charge Zi 
at rest in the EG at origin. The density induced in a homogeneous EG by a wave- 
vector and frequency dependent external perturbation ipexti^,^) plus the static 
potential of the charge is 



oo ^ 

n—lXn 

{q,qi, ...g„_i) 



X4>extiq - qi)--(t>ext{qn-2 - qn-l)(t>ext{qn-l), (A. 7) 

(l)ext{q) = ll^extiq) + Vq6{uj). (A. 8) 

Hence, to express the linear response function of the combined system of the EG 
with the charge Zi at the origin, we must collect in Eq. (jA.7[) terms linear in ipextiq)- 
With the use of symmetries (jA.ip . we readily arrive at the result 

X(q,k,w) = xi(q,^^)'5(q-k) + 2w,_feX2(q,t^,k,w) 

oo „ 

w,k,w) 

X Vq^q^...Vq^_^-kdqi...dqn^2, 

where the dependence on the wave- vectors and frequencies has been written explic- 
itly. Hence, by virtue of Eq. (|A.6p 

(qo • v) (q„_i ■ v) 



V 



n=l 

d 

X — fn{qo,i^, •■•,qn-i, w)dqo...(iq„-i 



j=0 



or 

oo n— 1 



(qo • v) (q„_i • v) 



It J. /• / 

P=(2.flm5:nE/^ 

n=l j=0 

X /ni(qo, ■•■,qr^-l)rfqo•••rfqn-l■ 
From Eqs. (jA.ip and (|A.3P it follows that /„ functions satisfy the same symmetry 
relations as x„ 



/„(go, ■■■qj,q3-qn-i, ...,qj-qj+i),0 <j< n-1. 



(A.9) 
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A.l. 



For brevity, throughout the derivations below by J... we imply the integral 
/ dqi...dq„_i.... By Eq. (jA.9p . for j < i we can write 

y qj/m(qo, •■■,qn-i) = - J <l]fn,n-i{q.o, ...qj-i,qj-i - q„-i, ...,qj-i - qj) = 

-J qj7ri,n-i(qo, ...qj-i,q„-i, ...,qj+i,qj) = - ^ q^ /„,«-! (qo, ...,q„-i), (A.IO) 

where we have used the integration variables q^, qn-i substitutions leaving other 
variables intact. But for j < n — 1 

J qi/n,Ti-i(qo, qn-i) = — J qj/ji,n-i(qo, ••■,qn-2,qn-2 — qn-i) 

= - J qj/n,n-i(qo, ■••,q«-i), (A.ii) 

which shows that for j < n — 1 

J qj/n,n-i(qo,...,qn-i) = 0, (A.12) 

and then by Eq. (fOO)) 

y qj/m(qo, ■••,qn-i) = 0, j<i. (A.ia) 

Further, if < j < i < n — 1, we can write 

J qi/nj(qo, •■■qfc, qfc+i, qn-i) = -J qi/n,n+fc-i(qo, ■•■q/t, q/t — qn-i, qfc - qfc+i) 

/(q.-q.O/,.n+.-,(qo,...q.,q„-l,...,q.+l)=/(q„-..-.-q.)/„, 
The first term in the last expression disappears due to Eq. (|A.13[) . and we have 

y qi/ni(qo, ■•■,qn-l) = J q.„ + k-tfn,n+k-]{(lo, ■■■,<ln-l)- (A.14) 

If i > j, then n + k — i < n + k — j, and using Eq. (jA.lSP again we have 

yqj/m(qo,---,qn-i)dqo---dq„-i=o, i ^ o,j ^ i. (A.is) 

Putting i — j and k = j — 1 in Eq. (|A.14[) . we obtain 

yqj/nj(qO,---,qn-l) ==yqn-l/n,n-l(qO,.--,qn-l), i^O. (A.16) 

A.2. 



I 



I 
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For < j < 7Z — 1 we can write 

/n — 1 ^ n— 1 

^(Is fn(){q.o,...,qn-i) = / ^(Is fno{q.o,.-,qj-i,(lj^i - q„_i, ...,qj_i -q^) = 

/ri— 1 „ n— 1 

[{n-j)cij-i - ^ q^] /„o(qo, •■•,qi-i,qn-i, ■•■,qi) = / [(n- j)qj_i - ^qs]/no(qo,•••,qrl-l)■ 
s=j s=j 

Then 

/n — 1 
X! /no(qo, q«-i) dqo.-.dqn-i = 

(n-j) y"qj_i/„o(qo,---,qn-i)dqo---cfqri-i- (A.17) 
Using Eq. (|A.17|) . we prove that 

y qj /no(qo, •■•,qn-i) = 

{n~ j) J C[r,-1 fno{<lO,---,q.n-l), j>0 (A.18) 

by induction from j = n — 1 to j = 1. For < i < n — 1, we can write 

n— 1 ^ n— 1 

q</«o(qo, •■■,qn-i) = / q» X! /"j(qo,qo - q„-i...,qo - qi) = / (qo - q*) ^ /ry(qo,qn-i---,qi) 

i=o i=o 

y" (qo - q«-j) X! /»j(qo,qi---,q«-i) = y"(qo - qn-j)/no(qo,qi---,qn-i) - y"q„-j/„,„-,(qo,qi...,q« 



J=0 

and then for < i < n — 1 



y"(qj + qn-j ~qo)/no(qo,---,qn-i) = -y q„-j/„,„-j(qo, ■•■,qn-i)- (A.i9) 

With use of Eqs. (fOS)) and (fOBl) . Eq. (fOO]) gives 

y"(nq„_i - qo)/«o(qo,---,qn-i) = - ^ qn-l/n,n-l(qo, qn-i)- (A.20) 



A.3. 

We write 



(27r)^ /■ 

P -Q ^ — 2~I'^X! X! / '^^ [("qn-1 - q,) • v] /„j(qo, ...,q„-i)dqo...dq„-i = 

^ n=l j=0 

,„ x3 cxD I n — 1 I 

— ^Im^y(qo-v)<|[(nq„_i-qo) • v] /„o(qo, qn-i)+X! [('^q«-i-qi) • v] /„j(qo, ...,q„_i) >dqo...dq„-i 



n=r I j=l 
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(27r 



^Im^ / (qo-v){[(nq„-i-qo) • v] /„o(qo, q«-i) + (qn-i • v)/„^„_i(qo, q„_i)} dqo...dq„-i = 0, 

n— 1 



where the second equahty is due to Eqs. (|A.15|) and (|A.16p and the fourth equahty 
is due to Eq. (|X20)) . 

Appendix B. Proof of Eqs. dU - dU) 

We can rewrite Eq. pop in the form 



Q = 



— ,0) — 



V V ■ ' duo 

or with use of the static sum rule (jC.ip of the next section 



w=0 



X(r"',r',0)drdr'dr"dr'", 



Q = 



Taking use of Eq. (fT5|) , we have 



[Vrno(r) • v] [Vr'no{r') • v] 9Imx 



dr dr' . 



or 

Q 



[Vrno(r) -v] [Vr'no(r') • v] 



j=0 



dr dr' , 



[Vrno(?') • v] [Vr'no(r') • v] 9Im/j;c(r, r', cj) 



9cj 



dr dr'~ 



j=0 



[V,no(r).v] [V,mo(r0.v] , 9 Imxifs(r", r'", 



x:^5(r"',r',0)drdr' dr" dr'" 



L)=0 



FinaUy, using the static sum rule (|C.2|) to simplify the second term, we immediately 
arrive at Eqs. © - (gl). 

Appendix C. Static sum rules for the scalar response-functions 

A static shift of an ion by an infinitesimal vector Ar must result in the same shift 
of the ground-state electron density. The perturbation corresponding to this shift is 

Ayo(r) - Ar- VFo(r) 

and the change in the electron particle-density due to this shift is 

A7io(r) = Ar- Vno(r), 

which leads us to the static sum-rule 

X(r, r', 0) V'yo(r') dr' = Vno(r). (C.l) 

The same sum rule evidently holds for the non-interacting KS density-response 
function 



XKs{r,r',0)VVKs{r')dr' = Vno(r). 



(C.2) 



I 



I 
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Finally, inverting the above two equations and using Eq. (jlSp . we arrive at the 
static sum-rule for the scalar xc kernel 

/,e(r,r',0)V'no(r')rfr' = VV^r). 

We note, that although the sum rules of this Appendix follow as the static limit 
from the corresponding dynamical sum rules of Ref. , the former hold for extended 
systems as well, while the latter are valid for bounded (finite) systems only, fact of 
which necessitates the independent derivation of this Appendix. 

Appendix D. Proof of the equivalence between Eqs. ([3]) and {[T]) 

The imaginary part of the KS response- function can be written as 

ImxK5(q,k,c.) = j [/(e,) - /(e,)](s-| e^^"" |p+) 

x(p+| e-*"^-- |s-)5(cc; - + e,) dpds, (D.l) 

where I p" *"), |s^), and Cp are the single-particles outgoing scattering states, incoming 
states'^, and the energies in the potential Vks{t), respectively, and / is the Fermi 
function. Equation (jD.l|) can be rewritten as 

IniXKs(q,k,u;) = "(^^ j " f{^s + i^)] x 
(s"| e^i"- |p+)(p+| e-^'^'- \s-)d{uj -ep + e,)dpds, 
and expanded at small lu to 

ImXi^s(q,k,w) = J ^i<^P ~ <^f)S{€s - ep) 

x(s-| 6*1'' |p+)(p+| e-*"" |s-) dpds. (D.2) 
Now we can write by virtue of Eq. ([3]) 

^^j5{ep~eF)Sies~eF){s-\[vVVKs{r)]\p+) 

x{p+\[w-\/VKsir)]\s~)dpds, 

where the square brackets in [v • VVR-5(r)] denote that the gradient applies to 
VKsi'<^) only, leaving the wave- functions intact. We can further write 

(s-| [v • VVKsir)] |p+> = (s-| V yVKsir) |p+> - (s-|y^f5(r)v • V|p+) = 

(s-f Go^Kss^lv VyKs(r)|p+) - (s-| VVs(r) v V |p + G+1/ksP+), (D.3) 
where we have used the Lippmann-Schwinger equations 

P^)=p) + G^VksP^)- 
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From Eq. (|D.3p we have 

(s-| [v • VVKsir)] |p+> = • s (s| Vks |p+) - i v • p (s" | F^^s |p), 
and using formulas for the transition-matrix ISU 

t(s ^ p) = (s| Vks |p+) = {s-\Vks |p), 

we have 

(s-| [v • VVKsir)] |p+) - i V • (s - p) i(s ^ p), 

and similarly 

(p+l [v • VVKsir)] \s-) = V (s - p) <*(s ^ p). 
We can then write for the noninteracting-electrons part of the friction coefficient 

Qi^^J 6iep~eF)S{es~eF)[w{s-p)]^\t{s^p)\^dpds. (D.4) 

Recal ling the expression of the differential cross-section through the T-matrix ele- 
ment ^1] 

da{dsp) 



(2^rii(s^p) 



and performing some integrations in Eq. (|D.4p explicitly, we arrive at Eq. ([T]). 

Appendix E. A proof that for a finite system tlie tensorial sum rules of 
Eqs. ([30 |) -([32 l) lead to the scalar sum rule of Eq. ([TS]) 

Equation ()32ll can be rewritten as 

c 

The next step, which involves integration by parts, requires the response function 
to vanish at infinity and, therefore, it applies to bounded systems only. For the 
latter case we can write multiplying Eq. (jE.ip scalarly from the left by V and using 
Eq. (HH) 



X»fc(r,r',c^)V'feV'Vo(r')dr' = c / x.fe(r, r', c^) V'^r' dr' - no(r) (E.l) 



X(r,r',c^) [v;-Fo(r') - lu^'^] dr' = Vjno(r), 

and after inverting 

X-\v,v',Lo)V,noir')dv' = V, Vb(r) - cc-V,. (E.2) 

A similar relation holds for xks 

J x^5(r,r',c^)V,no(r')dr' = VjVKs{r) - c^V,. (E.3) 

Subtracting Eq. (|E.2p from Eq. (|E.3[) and using the definition of Eq. we im- 

mediately arrive at Eq. [compare with Ref. . 
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Appendix F. Explicit expressions for the KS current-density response 
function, the xc kernel, and the x vector of Eq. {\22\i for a system with 
spherical symmetry 



F.l. KS current- density response function 

The KS response function can be explicitly written as 

E [C(r)V.#(r) - ^^(r)V.C(r)] 

X [rp{r')V^Mr') - ^o(r')V>^(r')] , (F.l) 

where "^oSj) and are the single-particle wave- function and eigenenergy, respec- 
tively, in the state a, and is the occupation number of this state. From the 
spherical symmetry of the problem it is easy to conclude that both xks and j^c 
leave invariant the subspace of the vectors of the form 

a(r) v + h{r) (r . v) r, 

where a{r) and h{r) are arbitrary scalar functions of r = |r|, and, hence, both x and 
the solution y of Eq. ([24|) are the vectors from the same subspace. Furthermore, by 
the definition p3p . y is the transverse vector, which property imposes the fulfillment 
of the relation 

o!(r) -I- 4 r h(r) -|- h' (r) = 0, 

and, therefore, Eq. (j24p effectively becomes an equation with respect to one unknown 
scalar function of the radial coordinate. With use of Eq. (jF.ip . we arrive at the 
equalities which are sufficient to evaluate xifS-dependent quantities in Eq. (j24p 



where 



c y X-ft:s.ii(i", I"', ^^)a(r', Lj)(5jfcdr' = a(r, tj)(5ifc + &(r,a;)rirfc, (F.2) 

c y" XKS,^j(r,r',w)V;-V;,F(r')rfr' = fc(r, tj) 5,^ + m(r, w)r,rfc, (F.3) 

d(r) = no(r)a(r) +(7(r), 
6(r,a;) = -q{r)/r^ +p{r)/r'^, 



and 



r 
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oo r 

X / dr' [G'(„.„+..(r,r')-G;;,,„_Jr,r')] y.{r') j a{r")dr" 





47rr^ 



a{r') -- J a{r")dr- 

OO 

p'{r) = 2rq{r) - r'[air)no{r)]' + ^ E (^'-^'^ ' ^^^'h-^ V'^^ 

oo 



(F.4) 



X J dr' 





1 



a(r') - / a(r")dr" 



27r^^" 



G,„e„+a.(r,r')-G,;,,_Jr,r')_ 

OO 

x(<5,„_i,;^ la + Si„,i,-ih)ya{r)J dr'[Gi^,,„+^{r,r') + Gl^,^_^{r,r')\ yc,{r')Ja{r")dr" , (F.5) 



*TO(r)]' 



OO 

X 1 [Gi,,,„+„(r,r')+G,;,,^_Jr,r')]y,„,feJr')F'(r')dr', 



where Vair) are the solutions of the radial Schrodinger equation 



dr'^ 



2VKs(r)+2e„ 



ya{r) = 0, 



and 



G;,.(r,r')=E 



i//.A:,,('-) j//.A:,,(r') 

e + ZTj - e/3 



is the Green's function. 



F.2. Exchange and correlation kernel 

In the case of spherical symmetry, we can write 

/a;c,ii(r,r',w) [a{r',uj)Sjk + b{r' , oj)r'jr'^.] dr' = d{r,uj)5ik + b{r,uj)rirk, (F.7) 



(F.6) 
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and, with use of the equations of Sec. 13.11 we have 



5(r) 



{3 [Airb{r)e^c{r)nn{r)+i£^a{r)no{r) {a' {r) + b' {r)) + r uj fi'^^{r) {r bi{r) + a[{r))] 



3raj2 no{r) 

+3c^C.c(r) ['irbi{r)+a[{r)+r^b[{r)] +ujfi,,{r) [7 a[{r) + r (106i(r) + r6'i(r) + 3a'/(r))]} 



b{r 



{ 12 z r2 b{r) noir) e^,(r) + 3 irnoir) a' (r) e'^^{r) +3tr^ no{r) b'{r) e',,{r) 



3r3 no(r) 

+ 12 r2 6i (r) CL(0 + 6i (r) ^;,(r) - 3 iu C^r) a'l (0 - fi^r) a'l (r) + 3 r Cc(0 a'l (0 

+r w ryi,(r) a'^ (r) + 15 C.c(r) 6'i (r) + 23 Cc> 77,e(r) fel (r) +?>t^u CAt) b[ (r) + 4 i^'^^r) b[ (r) 

-3ze,,(r)no(r) [a'(r) - r {5 r b' (r) + a" (r) + rH" (r))] + 3r lu C,,{r) a'Hr) + r lu fi,,{r) a'({r) 

SC.c{r)+^Vxc{r)] 6'i'(r)}, (F.8) 



where ai(r) = a(r)/no(r) and = 6(r)/no(r). 

We are using the LDA to the tensorial xc kernel of Sec. 13. li with the low- 
frequency expansion of /^^ ^(n, cu) and /^^ rp{n, ui) up to the first order in uj. Using 
the equalities I^U 



we have from Eqs. (|^^ 

C,xc{n,uj) = 0, 

ITi' 

fixc{n,uj) ^ — 



=0 3 dw 



= 0, 



(F.9) 



For /^^ jn(n, 0), we take use of its expression through the shear modulus Hx 



For the local density such that 1 < < 5 we obtain ^xc{n) by the interpolation 
between the values of Ref. ^^H. For rs < 1, we use the high-density approximation 



nkp{n) 
IOtt 



For 



a/^c,r(".'^) 



0=0 



, we are using Eq. (15) of Ref. 



mi 
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F.3. The X vector of Eq. 

Using the definition of Eq. ((22)) together with Eq. (|33p we can write 



x(w) = (rto - »^o)v + ^XKs(w)Vxx5(w)V(nov), 

where we have exphcitly written the frequency dependence of the functions. Ex- 
panding to the first order in lo 



X/<-5(^)V(nov) = Xks(O)^Kv) + w 



V(nov) 



and using the static sum rule (jC.l|) . we can write 

x(cj) = (no - no)v + ^xa's(w)V(v • Wks) 



— Xi^s(c^)VxKyO) 



c 
uj'- 



du! 



(v • VVks)- 



(F.IO) 



It must be noted that an expansion in uj of XKsi'-^) apphed to a gradient of an 
w-independent function starts from the term and, therefore, the second and the 
third terms in the above expression have leading terms of lu^ and to^, respectively. To 
evaluate Eq. (jF.lOp to the first order in u we derive and use the following equalities: 



^Xifs(w)V(v • Wks) = a{r, + b{r, lu) (v • r) r. 



(F.ll) 



9h{r,uj) - 15g{r,uj) 
a{r,uj) = — no[r), 



b{r,uj) 



Abglr^Lu) — 15h{r,Lu) 
167rr^ ' 



(F.12) 



h{r,uj) 




'hoirY 




■h,{r)- 








+ UJ 





(F.13) 



hoir) 
9o{r) 



1=0 



Gi{kF,kF,r) 



cos[(5/_i(fcF) - SiikF)], 



(F.14) 



hiir) 
51 W 



R 



oo 

irkpr^^-^ 

1=0 



Hi{kF,kF,r) 
Gi{kF,kF,r) 



{cos[,5;_i(fcf ) - ,5;(A:F)][5,'_i(fcf ) + <5,'(A:f )]- 



^sin[.._,(M-MM]}+{(|^-^ 
X sin[(5/_i(fcF) - 



Hi{ka,ki3,r) 
Giika,k(3,r) 



(F.15) 
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l'^yi,k^{r)yi-i.kf,{r) 

I r [yi^,,k, {r)vlk^ (r) - Vi,k^ {r)vU^u, (Ol- (F.16) 

Equations (jP.lip - (|F.16P are enough to explicitly evaluate the second term 
and to apply the first operator from the right in the third term of Eq. (jF.lOp . 
We, however, did not find a way to explicitly apply XksIO) ii^ the third term and, 
therefore, we invert xks{Q) on the set of the basis functions of Eq. (f34|) and finally 
use Eqs. (|E3| and ([F^ . 



Gi{ka,kf3,r) 
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